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Abstract 

For a class of large- N multi-matrix models, we identify a group G that plays the 
same role as the group of loops on space-time does for Yang-Mills theory. G is the 
spectrum of a commutative shume-deconcatenation Hopf algebra that we associate 
to correlations. G is the exponential of the free Lie algebra. The generating series 
of correlations is a function on G and satisfies quadratic equations in convolution. 
These factorized Schwinger-Dyson or loop equations involve a collection of Schwinger- 
Dyson operators, which are shown to be right-invariant vector fields on G, one for 
each linearly independent primitive of the Hopf algebra. A large class of formal 
matrix models satisfying these properties are identified, including as special cases, the 
zero momentum limits of the Gaussian, Chern-Simons and Yang-Mills field theories. 
Moreover, the Schwinger-Dyson operators of the continuum Yang- Mills action are 
shown to be right-invariant derivations of the shuffle-deconcatenation Hopf algebra 
generated by sources labeled by position and polarization. 
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1 Introduction 



Quantum Yang-Mills theory is at the heart of the microscopic description of strongly interacting 
particles. The group of based loops on space-time plays an important role in the formulation 
of Yang-Mills theory in terms of Wilson loops, which are gauge invariant variables containing 
much physical information pQ. Expectation values of Wilson loop observables are functions on 
this group. Expansion around the multi-color limit is a promising approximation method to 
solve Yang-Mills theory [2]. However, in the absence of a full-fledged differential geometry and 
analysis on the space of loops, progress in understanding and approximately solving the multi- 
color limit of Yang-Mills theory has been partly held up despite important early work j2[ [31 dj. 

On the other hand, hermitian multi-matrix models may be regarded as toy-models for Yang- 
Mills theory. The N x N matrices may be thought of as gauge fields at various space-time 
points, where N is the number of colors. It is natural to ask whether there is an analogue of 
the group of loops such that multi-matrix correlators are functions on this group. If so, can we 
interpret the factorized Schwinger-Dyson loop equations (fSDEfl of a matrix model in terms of 
the differential operators and products naturally associated to this group? Doing so may open 
up new perspectives and approximation methods for the large- N limit of multi-matrix models 
and Yang-Mills theory. 

1 The fSDE are quantum-corrected equations of motion for a matrix model in the multi-color (large- N) limit. 
They are analogous to the Makeenko-Migdal equations, 4 of large- N Yang-Mills theory. N = 3 in nature. 
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In this paper we show that there is indeed such a group G associated to a large class of 
multi-matrix models. We construct it indirectly as the group of characters of the commutative 
shufhe-deconcatenation Hopf algebrad- The latter is defined using the shuffle and concatenation 
products and reversal of order of matrices in a correlator. These are analogous to point-wise 
products of functions of loops, concatenation of loops, and reversal of loop orientation. In the 
simplest case of a single matrix and real- valued characters, G is the multiplicative group of non- 
zero real numbers. More generally, G is identified with the exponential of the free Lie algebra. 
We develop the rudiments of differential calculus on G using algebraic operations in the Hopf 
algebra. We find a large class of (formal) matrix models that can be formulated in terms of 
this group. We show that their fSDE are quadratic equations (in the convolution product) 
for a function on G. Moreover, there is one equation for each linearly independent primitive 
element of the Hopf algebra of functions on G. The Schwinger-Dyson (SD) operators, one for 
each linearly independent primitive, are shown to be right-invariant vector fields on G. Thus, 
given a prescription of which right-invariant vector field to associate to a given primitive, we 
can write down a system of fSDE for any group. For the group of relevance to matrix models, 
this prescription is encoded in the action. We find a large class of admissible actions and 
their SD operators, which include the gaussian, Chern-Simons and Yang-Mills matrix models as 
special cases. Finally, the SD operators of the Yang-Mills action are obtained and shown to be 
right-invariant derivations of the shuffle-deconcatenation Hopf algebra on a continuously infinite 
number of generators labeled by space-time position and polarization. However, in the case of 
Yang- Mills theory, we still need to pass to a quotient of this Hopf algebra to account for gauge 
invariance and recover the group of loops [H [6] (or need to gauge-fix and introduce additional 
generators for ghosts [7]), before we can look for physical solutions to the fSDE. 

In [8] the fSDE were formulated as conditions for the extremum of a large- iV 'classical' action 
(the Legendre transform of the entropy of operator- valued random variables). This viewpoint 
applied to generic multi-matrix models and also provided a variational approximation method. 
Here we develop a quite different group theoretic formulation, which only applies to a sub- 
class of matrix models. The distinction is very roughly analogous to that between the generic 
classical mechanical system and one whose configuration space is a group. However, due to this 
restriction, we find additional structures which closely mimic those present in Yang-Mills theory. 
We wanted to study these structures since they form the basis for an approximation scheme for 
multi-matrix models proposed in j9]. We hope our group theoretic formulation allows for a 
generalization to more familiar groups, where solutions to the fSDE may be more easily found. 
Our work continues the developments in the physics literature due to Migdal and Makeenko 
[UCEO], Polyakov [llj . Cvitanovic et. al. [12] , Gambini et. al. [HE], Tavares [6], Rajeev et. 
al. [HI [8], and builds on our previous papers [9J [7]. There are of course other approaches to 
multi-matrix models such as those related to integrable models and algebraic geometry, see for 
instance [IU [16] . 

By a A -matrix model |17} (9J [18] we mean a statistical system whose variables are a 
collection of random hermitian N x N matrices Ai , 1 < i < A with partition function Z = 
J dAe~ N tr . The integration is over all independent matrix elements. The action is the 
trace of a polynomial in the matrices tr S(A) = tr S 1 Aj where^l S 1 are the coupling 'tensors' 

technically, G can be thought of as an analogue of the group of generalized loops or extended loop group 
studied by Tavares[B] and Bartolo, Gambini and Griego [13]. We also identify a subgroup of G which can be 
regarded as an analogue of the smaller group of loops on space-time, see sec. |3] 

3 Capital letters denote multi-indices and repeated indices are summed. If I = iii2 • • ■ i n then S 1 = S' 1 and 
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(see sec. [5] for examples). The Ai model gauge fields at a collection of space-time points labeled 
by i and are N x N matrices in color space. Observables 3>/ = j^Aj are invariant under the 
global U(N) action A$ — > UAiU^. We are interested in their expectation values in the large- 
limit, the cyclic tensors (0 denotes the empty word) 

Gi satisfy a closed system of factorizecEl Schwinger-Dyson equations, conditions for the in- 
variance of Z under infinitesimal non-linear changes of integration variable. The latter are 
infinitesimal automorphisms of the tensor algebra generated by Ai , L\Aj = SjAj . The fSDE 
one for each letter i and word /, relate a change of action to a change in measure 

S Jlij2 G JlIj2 =8™*G h G h . (2) 

If we define A non-commuting sources £ l , we can form the generating series of correlators 
G(0 = EiGit 1 . Then! the fSDE can be written S*G(£) = G(g)?G(g). The SD operators 

S l = Y,(n + l)S^- 3 "D jn ..-D n (3) 

n>0 

are expressed in terms of left annihilation operators Dj which satisfy Dj^ n "' ln = Sj 1 ^ 12 '" 1 " or 
equivalently, [DjG]j = Gji. For more details on the fSDE, we refer to [8j[9l[18]. In this paper 
we do not have anything to say about the convergence of matrix integrals. We only use them 
as a formal device to generate the fSDE, whose structure we wish to investigate. 

2 Hopf algebra structure on correlations 

The space of based oriented loops 7 on space-time (up to equivalence under backtracking or 
retracing), plays a basic role in the Wilson loop formulation of Yang- Mills theory. This loop 
space forms an infinite dimensional non-abelian group, with successive traversal of loops as 
product 71 72 and reversal of orientation 7 as inverse. The information in this group of loops 
can be encoded in the algebra of (complex-valued) functions defined on it. Wilson loop functions 
W(j) = tr Pexp <f J 4 /J (7(t))7 /i (t)dt (trace of holonomy of the gauge connection A^{x) around 
the closed loop 7^(t)) form an adequate class of functions for this purpose pQ. Since the 
underlying loop space is a non-abelian group, the algebra of functions has the additional structure 
of a commutative but non-cocommutative Hopf algebra (under suitable hypotheses, this is a 
general property of the algebra of functions on any group [19 ). The point-wise product is 
(W1W2X7) = W\ (7)1^2(7) 1 the coproduct (AW) (71, 72) = W(7i72) encodes the concatenation 
of loops and the antipode (SW)(7) = W(7) encodes the inverse. The product and coproduct 
define compatible algebra and coalgebra structures, while the antipode turns this bialgebra into 
a Hopf algebra. Up to some technicalities, the underlying group of loops can be recovered as 
the spectrum (group of characters) of this Hopf algebra [191 125| [6]. 

For a multi-matrix model we did not know the analogue of the group of loops, but we 
did notice a bialgebra structure on the multi-matrix correlators Gi in connection with some 

Ai = Ai 1 Ai 2 ■ ■ ■ Ai n is a product of n matrices. The Kronecker Sj is equal to one if / = J and zero otherwise. 
4 Factorization[TO] is the property ("l , /i < l > / 2 • ■ • 4>i„) = Gi t ■ ■ ■ Gi n + O(jp) as N — » 00. 
5 Juxtaposition ^ Jl ^^ 12 denotes concatenation £ /l8i2 . 



3 



approximation schemes for the loop equations |9j. We recall this bialgebra structure and then 
define a compatible antipode to obtain the shuffle-deconcatenation Hopf algebra, which is the 
analogue of the Hopf algebra of Wilson loop functions. In the next section we will extract 
the underlying group from this Hopf algebra. The shuffle-deconcatenation Hopf algebra has 
appeared previously in other contexts [201 [T21 [B] . 

Let G(£) = J2i Gj^ 1 denote the generating series of multi-matrix correlators in the large- N 
limit. G(£) is an element of the vector space C((A)) of formal complex linear combinations of 
words £ J in generators chosen from the alphabet A = 1 < i < A} consisting of sources £*, 
one for each matrix Ai . The commutative shuffle product sh of two such series is (F oG)(^) = 
Y!,i(F ° G)j^ 1 where (F o G)j = J2i=juk FjGr ■ The sum is over all complementary order- 
preserving sub-strings J and K of I. For example, 

(F o G)ijk = F$Gijk + FiGjk + +FjGik + F^Gij + FijGk + F^Gj + Fj^Gi + F^G^. (4) 

Physically, the shuffle product is the product induced on gluon correlations by the point-wise 
product of Wilson loop expectation values ((WiW2)(7)) = (^1(7)) (^2(7)) > when path ordered 
exponentials are expanded in iterated integrals of gluon correlations in the large- N limit [9]. 
C({A)) with the shuffle product is the shuffle algebra on A generators Sh\. The empty word 
1 is a unit element for sh, with loF = Fol=F for all F £ Sh\ . 

Concatenation is defined as £ £ = £ , which extends linearly to (FG)(£) = FjGj£ ij . 
Using the inner product on C((A)) for which £ J form an orthonormal basis = 5 I,J , 

we can define the adjoint of concatenation or the deconcatenation coproduct, A = cone 1 by 
(F, GH) = (AF, G (8) H) . On monomials, A^ 1 = Sj K ^ J (8) £ K ■ It is extended linearly to series 
AF = J2jk Fjk£ J ®£, K ■ A is not co-commutative, it mimics the coproduct on Wilson functions 
coming from concatenation of loops (AW)(7i,72) = ^(7172) • We showed in [9] that A is a 
homomorphism of sh. The homomorphism of sh, e : Sh\ — * C which picks out the constant 
term, e(F) = F$, is a counit. Thus (sh,conc^ = A, l,e) is a bialgebra, the sh-deconc bialgebra. 

To turn this bialgebra into a Hopf algebra, we define an antipode S : C((A)) — ► C((A)), 
that mimics the antipode (SW)^) = W(7) on functions on loop space. On basis elements, let 

S(£, T ) = where ■ ■ ■ i n = in • • ■ ^h- (5) 

Extend it linearly to 5(F/£ 7 ) = E/t-l) 1 ' 1 -^ 7 • For example, S(? liaia ) = ai2h . This comes 
from comparing the expansions of W(j) and W(7) in terms of gluon correlations. For S to be 
an antipode it must satisfy several conditions which are usually summarized in a commutative 
diagram (see j21j). (1) It must be a homomorphism of the commutative shuffle algebra. On 
basis elements, this is the requirement S^ 1 o £ J ) = S^ 1 ) o S(!; J ) or equivalently, 

E (-i) 1 * 1 ^ = E (-i) 1 ^- (6) 

IUJ=K IUJ=L 

This is indeed true. Riffle shuffling two card packs (I, J) preserves the order of each. So reversing 
the order of the result of the shuffle (each summand on the lhs) is the same as reversing the 
order of each card pack (I, J) and then shuffling them together (each summand on the rhs). 
The minus signs just come along for the ride; S would be a homomorphism even without them. 

(2) The next two conditions S(l) = 1 and eS = e are obviously satisfied. 
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(3) The most interesting requirement for S to be an antipode is its compatibility with 
deconcatenation and shufflqj 

sh(S ® 1)A = sh(l <g> S)A = le or equivalently, 
^(-Ij'Voe^ = °£* = *0 for all J. (7) 

Putting I = i\ ■ ■ ■ i n these are the conditions 

n n 
p=0 p=0 

We have not found any nice proof of this, though we verified it explicitly for n < 3 and observed 
a pattern of cancelations for higher n which leads us to conjecture that it is an identity. Cartier 
[19j mentions that sh-deconc must form a Hopf algebra on general grounds, though we would 
still like an explicit proof of (JT]). The minus signs in the definition of the antipode are crucial 
for this compatibility condition to hold. In the sequel we will assume this condition is satisfied. 



3 Matrix model analogue of the group of loops 

The sh-deconc Hopf algebra we described is a commutative but non-cocommutative Hopf al- 
gebra, so it should be the algebra of functions on some non-abelian group. Which group is it? 
In the case of Yang-Mills theory, the corresponding group is that of based loops on space time. 
Remarkably, there seems to be an analogue of this group for hermitian A-matrix models. One 
might speculate that it is a group built from U(N) or a free group on A generators (since the 
concatenation algebra of correlations is the free associative algebra), but this is not the case. 
Rather, we will construct it as the group of complex valued characters (also known as the dual 
or spectrum) of the sh-deconc Hopf algebra. One might suspect that the analogue of loops are 
words in the generators of the shuffle algebra; but words do not form a group in a simple-minded 
way. Nevertheless, we will associate a (family of) group elements to each word and show that 
they form a subgroup of the spectrum. In another direction, using a result of Ree and Friedrichs 
[22\ 123] . we will identify the spectrum with the exponential of the free Lie algebra. 

Consider the set of real/complex- valued characters % °f the shuffle algebra, which are not 
identically zero. These are linear homomorphisms from the commutative shuffle algebra to the 
complex numbers. Suppose F,G £ Sh then 

X (F o G) = x(F)x(G) and for a, &eC, X (aF + bG) = a X {F) + b X (G). (9) 

It follows that X (l) = 1 for all characters X . We will define a group structure on this set and 
call it spec(STiA) or spec A . Suppose xiC 1 ) = X 1 ■ Then the complex numbers x ■> which we call 
the character coefficients, completely specify the character. For any F 6 Sh, X (F) = X (Fj^ 1 ) = 
Fjx 1 ■ The x form a dual space to the F 6 Sh, which justifies the upper and lower indices. We 
can also think of a character as a formal power series x = X^i • The identity is taken as the 
counit e : Sh — ► C defined as e(F) = Fq, which is a rather trivial homomorphism. In terms 
of coefficients, e 1 = Si. The product is non-abelian in general and is defined using A = conc^ . 

6 If C{{A}) were the algebra of functions on a group, these conditions would follow from the property that the 
product of a group element with its inverse in either order is the group identity. 
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So it encodes the monoid structure of concatenation. More precisely, yi/> = (x ® VO^, which 
is a map from Sh (8 S/i — > C (8) C . Then we identify C (8 C with C by multiplying the two 
components, to get a map Sh <8 Sh — > C. On basis elements, 

(xV0(£ J ) = x(e J ) iKe*) = <^xV- (io) 

It is extended linearly to the rest of Sh . So the product of characters is just the concatenation 
product of character power series. The formula for the product does not use the fact that 
characters are homomorphisms of shuffle. But we need the latter property to show that products 
of characters are also characters. A = conc^ is a homomorphism of the shuffle product as 
are x an d ip ■ Therefore, yV = (x ® VO^ * s a ^ so a homomorphism of the shuffle product: 
(xip)(F o G) = {xi>){F) (xVO(C)- Indeed, each side is equal to (we identify C (8) C with C) 

(AF)^ (AG) K>L flj* 4 UL X M (11) 

The inverse of a character is defined by composing with the antipode: = xS ■ S and x are 
homomorphisms of shuffle and so x 1 is a l so a homomorphism, and hence a character: 

x _1 (e J ) = xm 1 )) = (-i) |7| x(e f ) = c-i) 1 V- (12) 

The conditions XX~ 1 = X _1 X = e are precisely the same as the compatibility conditions (J7|) of 
the antipode S with product s/i and coproduct A . Indeed, using the homomorphism property 
of x an d the second equality above, 

(XX" 1 )^) = S!j K x^ J )(-l) lKl x(C R ) = x(#jc(-1) |J V t R ) = X(H) = (13) 

Similarly we verify that x~ l X = e - The shuffle algebra is the commutative algebra of functions 
with pointwise product on the group of characters. The value of a function F at the character 
X is obtained by evaluating the character on F: F(x) = x(F) = X 1 Fi ■ Moreover, (FoG)(x) = 
F( X )G(x) since x{F)x{G) = x(F o G) . 

We still need to find non-trivial characters. If x = X^i wants to be a character, x cannot 
be arbitrary. On the one hand, x{F) = Fix 1 may not converge, but we can consider polynomial 
F so that the series terminates. On the other hand, x must be a homomorphism of sh, and 
this imposes relations on the x 1 ■ F° r polynomial F and G, x(F ° G) = x{F)x{G) is satisfied 
iff xi^ 1 ° £ J ) = x(£ J )x(£ J ) f° r an 1} J or equivalently 

E X K = xV for all /, J. (14) 

/UJ=X 



These conditions were called the shuffle relations in another context [22] . They are the complete 
set of conditions for y to be a character. In detail, the first few shuffle relations are 

i 

X = 1, 

x ij + x ji = \V- 

yijk _|_ yJ*™ -|- = X*X 

y/i fc ' _l_ y^i' _)_ y» fe y _|_ y ki J l _)_ y^'- 7 + y feKj = x^X kl 

X ijkl + X m + X jkU + X jkU = xV'K e.t.c. (15) 

For rank n > 1 character coefficient tensors x n j there are [n/2] systems of linear shuffle 
relations (i.e. either ^(n — 1) or n/2 according as n is odd or even). The shuffle relations are 
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hierarchical, in the sense that the rank of the tensors on the lhs ( |/| + 1 J| ) always exceeds the rank 
of the tensors on the rhs ( \I\ and | J| ). So we can think of these as linear equations constraining 
the higher rank \ K m terms of the lower rank ones which appear quadratically as sources on 
the right. This structure is reminiscent of the matrix model fSDE: S Jllj2 Gj 1 u 2 = (5j ll/2 G/ 1 G/ 2 
Naively, we expect a large space of solutions to these constraints, since there seem to be 
a lot more degrees of freedom in the x 1 than there are shuffle relations. In particular, the x* 
are unconstrained. Regard x lJ as a matrix. Then its symmetric part is completely determined 
by the x% but its anti-symmetric part ^(x^ — X^) is not. 

In the case of loop space, we can find examples of characters easily. For example, given a 
loop 7(t)o<t<i on space-time M, its value on an element of F £ Sh(M) is j(F) = J F where 
the rhs is the iterated 'Chen' integral[5] |6] of the linear combination of tensor products of one 
forms. For example if F = a <S> for a pair of 1 -forms a and (3, then 

7 (F) = / dh f 1 dt 2 Qifrfa)) &(7(t2)) f(h) f{t 2 ). (16) 
Jo Jo 

For the shuffle algebra on a finite number of generators, we might imagine that the analogue of 
a loop is a word £/ and define a linear functional on Sh by £i(F) = Fj . However, this is not 
a character since = 5f , whereas for a character we must have x(l) = 1- Though single 

words are in the dual of Sh regarded as a vector space, they are not (with the exception of the 
empty word) in the dual of Sh regarded as an algebra. 

Thinking of a character as a formal power series x{0 = J2i X ! £i > we as k whether there are 
any characters aside from the identity e. To begin with, we show using the shuffle relations that 
there are no non-trivial polynomial characters. Suppose x ls a non-trivial polynomial character, 
of degree n — 1. What this means is that x K = for all words K of length \K\ > n > 2, 
but with x 1 ^ f° r some I of rank |/| = n — 1. Then consider the homomorphism condition 
xii 1 ° i 1 ) = X 1 X I which is the same as 

E X K = X I X I - (17) 

I,IUI=K 

The rhs is non-vanishing by assumption. But the lhs vanishes since it is a linear combination of 
character coefficients of rank 2n — 2 > n. Thus we have a contradiction. So the only polynomial 
character is the identity x = e - To find non-trivial characters, let us specialize first to the case 
of a single generator. 



3.1 Characters of the shuffle algebra on one generator Shi 

For A = 1 , a character is a formal series x = SnLo Xn£, n in one generator £ . The condition that 
it be a homomorphism of sh is xo = 1) an d the following shuffle relations for each Xn, n> 1: 

Xn = XrXn-r for r = 1, 2, 3, • • • , [re/2] . (18) 

In more detail, xo = 1, ^X2 = Xi, §X4 = xh ^Xa = X1X3, 5^5 = XiX4, 10x5 = X2X3 
etc. The general solution is a 1-parameter family Xn = ^\Xi fo r re > 0. We write x = 
In particular, there are no polynomial characters. Moreover, if xi 7^ V'l then e Xl ^ and e^ 1 ^ 
are distinct characters as they have different coefficients in their power series expansions. The 
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identity character is got by choosing %i = 0, in which case % = 1. The product x^P is the 
character whose value on monomials is 

W)(H = EXn-r^ (19) 
r=0 

(xVOi = Xi+V'i so xV' is the character e^ 1 " 1 "^ 1 ^ , which agrees with the usual rule for multiplying 
X = e Xl ^ and if) = e^ 1 ^ . The product is abelian, since we have a single generator. The inverse 
of x = eX1 ^ is X 1 = e~ X1 ^. We call the group of characters of Shi as spec(STii) or specx for 
short. Though space-time has been reduced, in a sense, to a single point, spe^ is a continuous 
abelian group parameterized by one real/complex number xi ■ Indeed we can even define a 
one dimensional abelian Lie algebra on the vector space {xi£ I Xi 6 C or R} with Lie bracket 
[xi£>Y'i£] = 0- X = eX1 ^ is an exponential map from the Lie algebra to the group. 

If we consider real-valued characters, then / : e x ^ i-> e x * is an isomorphism from spe^ to the 
multiplicative group of non-zero reals, R* . For complex- valued characters, / is a homomorphism 
from spec l onto C* , the multiplicative group of non-zero complex numbers. Its kernel is the 
subgroup generated by the character e 2nt ^ , i.e. the subgroup {e 27 ™" 1 ? | m £ Z} . 

We think of the shuffle algebra Sh\ as the commutative algebra of functions on the group 
specj^ . The value of F = F n ^ n at x is 

F{X) = X(E F nO = E XnF n or equivalently C(x) = x(C) = Xn- (20) 

n n 

The shuffle product is the same as the point- wise product of functions on spec 1 , since characters 
are homomorphisms of the shuffle algebra 

r°r)(x) = xr°r) = ( n+ n m \x(e +m ) = ( n+ n m )xn +m = xnx m = e(x)c m (x)- (21) 



3.2 Group of characters of Sh\: Pure characters 

We now discuss the group of characters of the sh-deconc Hopf algebra on A > 1 generators. 
Corresponding to the inclusions £j {£1, • • • , £a} we get A abelian one-parameter subgroups 
of spec(Sh\) for free, namely e x ^, e x • • • , e x ^ A . For instance, the value of e x ^ 3 on a 
basis-element of the shuffle algebra is 

( e x 3 C3)(tl\ = f if I = 333. ..33 (n times), and ^ 

V AS ' \ otherwise. V ; 

We call a character x pure if it does not mix the letters {£i,- • ■ , £a}, i-e. X 1 = whenever 
£ J contains at least two distinct letters. The only pure characters are the ones listed above. 
Another way to look at it is that given a letter ^ from the alphabet {£i, • • • , £a} , we get a 
1 -parameter family of pure characters e x &i . 



3.3 Mixed characters 



A natural question is whether given a word £/, it is possible to obtain a (family of) character(s) 
associated to it? This would in a sense be analogous to associating the character of the shuffle 
algebra Sh{M) , J F to each based loop 7 on space-time where F £ Sh(M) . 
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We call a character x mixed if it is not pure, i.e. if there is a word £/ containing at least 
two distinct letters such that x 7^ 0- A large class of mixed characters can be obtained by 
multiplying pure characters. Given a word £/ = ■ • • & n and a sequence of complex numbers 
{x 1 ) " " ' >X n }> we define a character via the product of pure characters (products of characters 
are characters (sec. S}) 

^ = gx 1 ^! e x% 2 . . . e x n &„ . (23) 

A word-sequence pair x) = (&i ■ • •&„) {x\ " • • X™}) determines a character. The inverse of 
X is also of the same form, = e - * 1 ^™ e~ x ■ ■ • e* n £n and satisfies t/^ -1 = ip -1 ^ = e. 

Reduced form of word-sequence pairs: Given a word Cii " " " Cij-iCiiCij+i ' " ' £i n anc ^ a sequence 
(x j • • • , X" 7- 1 ) X 3 iX^ + > ■ " ) X n ) ) we define the reduced form of the pair. If a pair of adjacent 
letters coincide, £^ = & . +1 , then delete & . and x J > an d replace x-^ 1 with + X"^ 1 to get 
a new word £ h ■ ■ • ^ J -_ 1 Ci j+1 • • • an d a new sequence x 1 , • • • , X J + X J+ \ • • • , X™- The 
resulting character is the same as the original one. Moreover, if any of the numbers Xk vanishes, 
just delete it and the letter £j fc . Proceeding this way, we get a word whose adjacent letters are 
always distinct, and a sequence of non-zero complex numbers (the one exception is if the word is 
empty). Such a word along with its sequence of complex numbers is in reduced form. Of course, 
the length of the word is the same as the length of the sequence. Two pairs (£/, x) an d (£j> VO 
are equivalent if they have the same reduced forms. Equivalent pairs correspond to the same 
character. For example, the reduced form of the pair (£3£i£i£2j {i, —3, n, 0}) is (£3^1 , {i, it — 3}) 
and corresponds to the character e l ^ 3 e^ 7T ~ 3 ^ 1 whose inverse is e^ 3 ~ n ^ 1 e~ l ^ s . The reduced form of 
the identity character is the pair consisting of the empty word and the empty sequence, (0, {}) . 
A pure character e x "1 has reduced form (£j 1 , {x 1 })- 

We can multiply two characters corresponding to the words ' " £i„ an d £,i n+1 • • • £« n+m 
to get a character corresponding to the concatenated word Ch' ' ' d„+ m an d the concatenated 
sequence of complex numbers x , " " " > x n+m '■ 



gX 1 ?*! gX 2 ?i 2 . . . gX n Cin gX n + 1 ?j n+1 gX n + 2 fi n + 2 . . . gX n+m £i n + m _ (24) 



This product is as non-abelian as it can be. The products of pure characters form a subgroup 
of the group of characters of Sh\ . From our construction, this subgroup is the free product of 
A copies of spec 1 

Group of prod, of pure charac. = spe^ * spe^ * • • • * spec 1 (A factors). (25) 

We know that the free group on n generators F n is the same as the free product of n copies of 
the integers F\ . By contrast, we will see that specx * specx * • • • * specx is a proper subgroup 
of spec^, at least if we restrict to finite products. It is interesting to know the appropriate 
topology for such free products of continuous groups (spec 1 = R* for instance). 

The exponential of any finite linear combination of generators x = e x ^ j is a character. To 
show this, consider 



X = e 



x 3 i 



oo -. 00 1 

E -^Hn ■ ■ ■ x J %„ = E ■ ■ ■ * in ^ • • • ( 26 ) 



n=0 n=0 



The coefficients are symmetric tensors x J ' 1 '" J " B = ^X"' 1 " " " X Jn which must satisfy the shuffle 
relations X IUJ = X I X J ■ Taking I = i\ ■ ■ ■ i n and J = ji ■ ■ ■ j m , on the lhs there are ( n ~^ m ) 
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terms all of which are equal, so 



LHS 



( 



n + m 



n 



) 



X 



^1 ' "^njl ' ' 'Jm — 



n!m! 



1 



x 11 ■■■x tn x 




RHS. 



(27) 



e x & is thus a character, but since £j do not commute it cannot be written as a finite product of 
pure characters. The inverse of x is X -1 = e~ xJ ^ , which is also of the same form. Since products 
of characters satisfy the shuffle relations, finite products of the form x = e^'^'i e x ' 2 ^2 • • • e* 1 ™^™ 
form a group which properly contains spe^ * • • • * spec A and is a proper subgroup of spec A . As 
before, we can put any such product in a reduced form. 

Exponentials of arbitrary non-linear polynomials in the generators are not characters in 
general. For instance, it is easy to see that x = e^ 1+ ^ 2 does not satisfy the shuffle relation 



X 12 +X 21 = X*X 2 • O n the other hand, the Baker-Campbell-Hausdorff formula (BCH) tells us that 
e x% e ¥ti = e x*ei-h^'<j+3xVKi,€j]+^xW*Ki.fe^k]]+- . Using BCH we can reexpress products of 



exponentials occurring in the above subgroup of spec A as exponentials of linear combinations 
of nested commutators of the generators . Aside from x l £i j these will be certain infinite 
linear combinations since the £j do not commute. This suggested to us that exponentials of 
finite (or other infinite) linear combinations of nested commutators of may also be characters. 
While we verified this in some simple cases, the calculations rapidly get laborious. So we were 
pleasantly surprised to find this proven in the work of Ree [22] using a theorem of Friedrichs 
|23j . More precisely, linear combinations of iterated commutators of the generators £j are 
called Lie elements. They are obtained using the operations of taking Lie brackets [., .] and 
linear combinations but not products such as Ci^j ■ F° r example. C*£j + C ljfc [£j, + 
C l ^ kl [[ii, [Cfc, 6]] are Lie elements for any tensors C\ C iJ , (7 U . Ree proves that exponentials 
of Lie elements are the only formal series satisfying the shuffle relations. In other words, the 
group spec^ consists precisely of exponentials of Lie elements. This characterization will be 
useful since the Schwinger-Dyson operators of an interesting class of matrix models related to 
Yang-Mills theory, turn out to be Lie elements. Now that we have identified the group spec(Sh\) 
which plays the role of the group of loops, we can formulate differential calculus on it. 

Technically, the group of loops on space-time is a proper subgroup of a larger group of 
generalized loops (in the language of [6] ) or extended loop group (in the language of [U H3] ) . 
spec^ is the analogue of this larger group, while the group generated by pure characters is the 
analogue of the smaller group of loops. Both in the space-time and matrix model settings, the 
larger group behaves akin to a classical Lie group and appears to be the correct physical setting 
for the fSDE. For instance, the SD operators of several matrix models can be interpreted as 
right-invariant vector fields on the larger group, but not on the smaller one (sec. [5]). Moreover, 
just as the group generated by pure characters is a free product of A copies of spec 1 , the group 
of loops is also a free product (free group generated by the based loops) . 

4 Differential calculus on the group G = spec A 
4.1 Functions 

By a function on the group G = spec^ we will mean an element of the shuffle algebra F(£) = 
Fj^ 1 . Its value at the character x = X T £,i £ spec A is given by F(x) = x(F) = Fix 1 ■ The ring 
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of such functions is the commutative shuffle algebra Sh\ with point-wise product of F and G 
at x given by the shuffle product Y^i{F ° G)ix I ■ 

4.2 Vector fields and the Lie algebra of spec A 

By a vector field V on spec^ we will mean a derivation of the shuffle algebra, i.e. a map that 
takes functions to functions V : Sh — * Sh that is linear over the complex numbers and satisfies 
the Leibnitz rule V{F o G) = VF o G + F o VG . This extends the concept of vector fields to 
settings more general than differentiable manifolds [24J . The derivations must form a left-module 
over the shuffle algebra as well as a Lie algebra over the ring of functions (i.e. FV must be a 
vector field and the structure functions of the Lie algebra of vector fields must be in Sh). 

From sec. 13.31 a necessary and sufficient condition for x to be in G = spec^ is that 
logx(C) must be a Lie element, i.e. a linear combination of iterated commutators of £1, • • ■ ,£a- 
The set of Lie elements is closed under commutators and forms the free Lie algebra of rank A 
(FLAa). It is a free Lie algebra since there are no relations besides linearity, antisymmetry 
and Jacobi identity that are satisfied by the commutator brackets. This is reminiscent of the 
exponential map from the Lie algebra to a Lie group. So we expect the Lie algebra of right (or 
left) invariant derivations of Sh\ (which should play the role of Lie algebra of the group spec A ) 
to be isomorphic to FLAa • We will see that this is indeed the case. 

In [9] we showed that linear combinations of iterated commutators of left annihilation A 
are derivations of the shuffle algebra, so they may be considered vector fields on spec^- We 
recall why A satisfies the Leibnitz rule: [Di(F o G)]j = [F o G]u = J2iiul 2 =il ^h^h ■ Now 
either i G I\ or i G h , so 

[A(^°G9]/ = E {F ih G h +F h G ih )= £ {[DiF] h G h + F h [DiG\ h ) 
hul 2 =l hul 2 =l 
= [{D i F)oG\ I ^[Fo{D i G)] I . (28) 

Moreover, commutators of left annihilation do not satisfy any relations besides linearity, anti- 
symmetry and the Jacobi identity. So iterated commutators of A span a FLA of rank A . There 
is a standard basis An for the FLA, labeled by Lyndon words L . A^L) is a particular iterated 
commutator of A; 's, where L = l\ ■ ■ ■ l n is a Lyndon word (see appendix lA.lh . So an element 
of the above FLA of vector fields is written V = V^Dr^ where Vl 1 are (real or complex) 
constants, and the sum is over all Lyndon words L . The value of such a vector field at the point 
X = X T ii is the 'tangent vector' V^x^^iL) = VaDn^, since x® = 1 f° r a character. These 
'constant coefficient' vector fields can be regarded as forming a sub-algebra of the Lie algebra 
of spec A - For, by evaluating at the identity x = e > they span a space of 'tangent vectors' at 
the identity. We believe these constant coefficient vector fields should be regarded as the whole 
of the Lie algebra of G. The structure constants of the FLA of basis vector fields An are 
denoted [Al)> Am)] = c l m^(n) > see appendix lA.il for examples. 

A vector field V on a group is distinguished if it commutes with the action of the group on 
itself by multiplication (encoded in the coproduct A on the Hopf algebra of functions). Roughly, 
this means A V must equal VA . But that cannot be quite right since A : Sh — > Sh <8> Sh while 
V : Sh — > Sh, so we must specify whether V acts on the first or second slot. In fact we must 
distinguish the right action from the left action, which leads to the definitions AV = (V ® 1)A 
and AV = (1®V)A |25| . Alternatively, let R g be the right translation by g on a group G and 
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R* and R g ^ the pull-back and push- forward maps. Then the push- forward of a vector field V 
acts on a function / according to Rg^Vf = R*~ 1 VR*f . Now a vector field is right invariant 
if R 9i V = V or in other words, R*Vf = VR* g f for all /. But R* is the pull-back induced by 
right multiplication in the group, and multiplication in the group is encoded in the coproduct A 
in the Hopf algebra of functions on G. This justifies the definitions of right and left invariant 
derivations by AV = (V (g> 1)A and AV = (1 <g> V)A. 

If vector fields V and W are right invariant, then so is their product VW (though not a 
vector field) and hence also their commutator [V, W] (which is a vector field, since commutators 
of derivations are also derivations). To see this, 

(V®l)A = AV (V ®1)AW = AVW. (29) 

Using the right invariance of W to re-express AW , this becomes 

(V <8> 1)(W <S> 1)A = AVW (VW (g> 1)A = AVW. (30) 

Thus right invariant derivations form a Lie algebra, which will serve as a substitute for the Lie 
algebra of right-invariant vector fields on the group G . 

It remains to identify the right-invariant derivations in the case G = spec(Sh\) . We will 
show that the constant coefficient vector fields V = V^Dns are in fact right invariant. It is 
straightforward to check that Di is a right-invariant derivation: 

(A®1)AG = (Di®i)Gj K S J ®£ K = GijKS J ®t K 

A(AG) = A{G iI i I ) = G iI 5 I JK e ®t K = G lJK i J ®i K - (31) 

By the previous result, we deduce that iterated commutators of Di are also right-invariant. 
For example, [[Di, Dj], [D^, Di]] and [Di, [Dj, D{\] are right invariant derivations. Moreover, 
the condition AV = (V (g> 1)A is linear in V , so real/complex linear combinations of iterated 
commutators of Di are also right-invariant. In other words, constant coefficient vector fields 
V = V*D(r\ are right-invariant on G. In particular, the Schwinger-Dyson operators S % = 
Ag lk g jl [Dj, [D k ,D{\] of Yang-Mills matrix models with action S = tr g lk g jl [Ai , Aj ] [j4^ , A[ ] are 
right invariant vector fields (see sec. [5]). 

A more general derivation of the shuffle algebra is obtained by allowing non-constant co- 
efficients (elements of Sh), V = V L (£)Dr£ l \ = Vj- 1 ^ 1 Di L \. V acts on a function as VF = 
V L (0 (D(l)F{0) • The derivation property follows from that of -D(l) and commutativity of 
the shuffle product 

y L (0 D {L)(F oG) = V L {£) (d {l) F oG + Fo D {l) g\ = VF o G + F o VG. (32) 

Thus, we can think of V = V L (^)D^ as a vector field on the group G = spec A . It is also clear 
that these derivations of Sh form a left module over the ring of functions on spec^ • Indeed, 
shuffle multiplying a derivation V L (^)D^ from the left by a function F(£) gives another 
derivation {F(£) o V L {t)}D (L) . 

It is interesting to know whether there are any more derivations of the shuffle algebra; these 
seem to be adequate for us. We call the space of variable coefficient derivations VectiG) . In 
a sense, the -D(l) form a globally defined moving-frame so that the tangent bundle of G is 
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trivial, i.e., spec^ is parallelizable just like any Lie group. Vect(G) forms a Lie algebra with 
Lie bracket (all products are shuffle products) 

[V, W] = (VW L )D {L) - (WV L )D (L) + cl M V L o W M D {N) . (33) 

Here VW L is the action of the vector field V on the function W L (£), VW L = V M (£) o 
D(m)W l (£,) . We used the fact that Drn satisfies the Leibnitz rule with respect to sh. This is 
the analogue of the formula for Lie brackets of vector fields on a manifold [v,w] = v l diiu J dj — 
w l diV^dj. Dr L \ play the role of di, except they don't commute. 

The value of vector field V = V L {^)D^ at x £ G is obtained by evaluating the coefficient 
functions on the character x> to get a 'tangent vector' Vfx 1 ^(L) ■ I n particular, the value of 
a vector field at the group identity \ l = &\ is the 'tangent vector' V^D^ . Thus the space 
of tangent vectors at the identity is the same as the space of constant coefficient vector fields, 
which we have also observed to be right-invariant vector fields. Their Lie algebra is isomorphic 
to FLAa ; we are now justified to think of it as the Lie algebra of the group G = spec A . 

We argue that non-constant coefficient vector fields V cannot be right-invariant. Suppose 
there were right-invariant V , with Vj L for some non-empty I. V as well as the distinct 
right-invariant vector field with constant coefficients, V^Dus , both evaluate to the same tangent 
vector at the identity: V^D^ . However, there should be a unique right-invariant vector field on 
G that is obtained by right translating the tangent vector Vl'Dij^ , and the constant coefficient 
vector field V^D^ serves that purpose, so V could not be right-invariant. We would still like 
a combinatorial proof that non-constant coefficient vector fields cannot satisfy AV = (V® 1)A. 

4.3 One-forms 

A 1-form is a linear function from Vect(G) to Sh. The dual basis of 1 -forms 9 L is also labeled 
by Lyndon words L. On basis vector fields 6 [D<m)) = an d extended linearly to Vect(G), 
QL^yM (QD(\fi) = V L {^). A general 1-form u> = lol(£)9 l is a linear combination of the basis 
9 L with coefficients coming from Sh. Though Di L \ was defined (app. [AJ through iterated 
commutators of left annihilation D{ , we have not built 6 from 9 l . The exterior derivative of 
an element of the shuffle algebra is a 1-form defined by its action on vector fields 

dF(V L (OD [L) ) = V L {0 o D {L) F. (34) 
From this we can read off that if dF = (dF)i9 L , then the components (dF)i, = D^F. 

4.4 Differential calculus on spec^ 

Let us illustrate the above formalism in the simplest case of one generator. The spectrum of 
the shuffle algebra on one generator consists of the exponential series x(0 = S n >0 Xn£, n with 
Xn = rrX? • Functions on specx are elements of the shuffle algebra F = F n ^ n with the value 
Fix) = Sn>0 FnXn ■ = is the only left annihilation operator and its commutators 

vanish. So the Lyndon basis of the free Lie algebra on one generator is just D and it is an 
abelian algebra. Moreover D is right invariant since 

AD£ n = A^- 1 = Z P ®Z q and 

p+q=n— 1 
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p+g=n p+g=n 



(35) 



are equal. The general vector field is V = V(£)D where V(£) = J2n>o Vn£ n and it is right 
invariant iff V(£) is a constant. V restricts to the tangent vector J2 n >o VnXnD at the point x 
on the group spec 1 and to the tangent vector VqD at the identity. 

The 1-form dual to the vector field D is 6 with 9(D) = 1. The general 1-form is uj = u(£)6. 
The value of uj on the vector field V(£)D is the shuffle product uj{V) = uj(£) ° V(£)- The 
exterior derivative of a function is dF = (DF) 9 and for a monomial, d^ n = £ i n ~ 1 9 . The 1-form 
uj restricts to the co- vector J2 n >o U) nXn9 at the point % an d to the co- vector ujq9 at the identity. 

5 Factorized Schwinger-Dyson equations on G = spec A 

We found the group spec A that plays the role of the group of based loops, and identified the 
rudiments of differential calculus on it. Now we'd like to formulate the factorized Schwinger- 
Dyson equations of matrix models in terms of spec^. This may indicate how to generalize 
the fSDE to groups we are more familiar with, and thereby provide more insight into their 
solutions. The fSDE 5*G(£) = G(^)^G(^) are a system of equations for the moment generating 
series = Gj^ 1 . G(£) is an element of the shuffle algebra or equivalently, a function on 

the group spec^. S l are called the Schwinger-Dyson operators. The only a priori conditions 
are that G must evaluate to 1 at the group identity x = e (i- e - G$ = 1), the coefficients Gj 
must be cyclically symmetric and satisfy the reality condition G} = Gj. These follow from the 
physical requirements of normalization of expectation values = 1, cyclicity of the trace 

Gj = lim^r^ 00 (|^^4/) and reality of the matrix model action S(A). It will be useful to keep the 
gaussian, Chern-Simons and Yang-Mills (g %3 = g 3 " 1 ) examples in mind, 

S G = I tr C^AiAj , S cs = tr C i3k Ai [Aj , A k ] , S YM = g lk g jl [A , A,-] [A k , A t ] . (36) 

Their SD operators were obtained in [9]: 

4 = CtWj, Sh s = >/=Tk(C^* - C*y)[D fcj £),-], S YM = 4g lk g 3l [D 3 , [D k , A]]- (37) 

There is one fSDE for each letter £ l . But what does a letter mean in terms of the group? We 
noticed that letters are primitive elements of the commutative Hopf algebra of functions on 
G . Primitive elements P are those that satisfy AP = 1 (g) P + P ® 1 , where A = conc^ is 
deconcatenation, which is not cocommutative^- On monomials, A^ 1 = 5j K £ J '®£ K ■ Let us show 
that the only primitive elements are linear combinations of letters £ l , 1 < i < A. For an element 
= Pj^ 1 of the shuffle algebra to be primitive, we need Pjk£, J <8> i K = Pii 1 <8> 1 + 1 <8> PiS, 1 ■ 
This is equivalent to the requirements P$ = and J2j,k^® Pjk£ J ® £ K = 0. These conditions 
are satisfied iff P$ = and Pj = for |/| > 2. In other words, linear combinations of letters are 
the only primitives. So if we pick any basis for the vector space of primitives (such as the letters 
themselves), we will have one fSDE for each basis element. This is again consistent with the 
fact that we could rewrite the fSDE as wlS % G(<£) = G(£)u>;?£*G(£) for any non-singular A x A 

7 In other words tA 7^ A, where t(F Cg> G) = G ® F reverses the order of factors. Note that this is distinct 
from the Poincare-Birkhoff-Witt construction where the primitives of a cocommutative Hopf algebra form a Lie 
algebra whose universal envelope is the Hopf algebra. 
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matrix Wj . This characterization in terms of the primitives of the Hopf algebra of functions on 
G applies to any group. 

The rhs of the fSDE G(£)£ l G(£) = GhGj^ HiJ involves concatenation. It can be under- 
stood in terms of the convolution product of functions on the group G. Given a non-abelian 
group G, there are two natural dual Hopf algebras associated to it, the commutative algebra 
C G of complex functions F(g) on it with pointwise product (FG)(g) = F(g)G(g), and the 
non-commutative group algebra CG = {X) 9 gG F(g)g} with convolution product (FG)(g) = 
Y^heG F(K)G(h g) . The coproduct of the first becomes the product in the second and vice versa 
using the duality {J2heG F(h)h, G) = Z) s gG F(g)G(g) . In our case, sh-deconc is the commuta- 
tive Hopf algebra whose coproduct is deconcatenation. The dual Hopf algebra cone- deshuf fie 
is the convolution algebra of functions on G, whose product is concatenation {FG)!^ 1 = 
FjGk£ JK ■ Thus is the convolution of with the primitive element con- 

volved again with This formulation again applies to any group. 

For each primitive element £ l , the lhs of the fSDE is the SD operator S l acting on the 
function G(£) on the group. In the Gaussian, CS and YM examples (|37|) , S l is a complex-linear 
combination of iterated commutators of left annihilation. Let us restrict attention to models 
where this is the case. Then from our discussion in sec. 14.21 we conclude that S 1 are right 
invariant vector fields on G . For this we think of the right-invariant vector fields (or Lie algebra 
of G ) as represented linearly on the space of functions on the group. 

So far, we have formulated practically everything in the fSDE in terms of concepts that 
generalize to any group without reference to matrix integrals. It only remains to specify which 
right-invariant vector field S l to associate to a given primitive element . For this we need 
additional data beyond the mere specification of a group G. At present we do not have a 
prescription of which right invariant vector field S l to associate to a given primitive £ l , that 
would apply to an arbitrary group. But we have a rough idea. The additional data is the 
specification of an action, and the prescription is the passage from action to SD operators. 
However, even for the group spec^ of relevance to matrix models, we only gave 3 examples 
(j37j) of actions leading to SD operators which are right-invariant vector fields for each primitive. 
Moreover, we know that many actions do not lead to SD operators that are right invariant vector 
fields (eg. S = tr A 4 ). So we postpone a general characterization of admissible actions and 
their passage to SD operators, which would apply to any group. Instead, we seek more examples 
of matrix model actions whose SD operators are right invariant vector fields on spec^ . 

Given a matrix model action S(A) = tr S 1 Aj , we define S(G) = S J Gj with cyclic S J and 
Gj . The SD operators S l are obtained from the variation in the action under the infinitesimal 
variations SAj = L\Aj = SjAj , i.e., by applying^ L\ to the actiorj§ 

S*G(0 = eL\(S J Gj) = eS J 5f ij2 G JlIj2 

= ^'S ''•'■'■'( = ^T l \Ji\S iJ GjiZ I = ]T | Ji\S iJ (DjG)ie. (38) 
J, i j,i 

From this we read off S l (which are not right-invariant vector fields in general) 

& = \Ji\S U Dj =J2(™ + l)S ij ^D jn ■ ■ ■ D n . (39) 

J n>0 

8 L\ are infinitesimal automorphisms of the tensor algebra [8] [9]. Their action on Gj is L\Gj = 8j 1% 2 Gj 1 u 2 . 

9 | J is the length of the word J and J is the reversed word. Cyclicity of S J and Gj are used in the fourth 
equality. Dj is left annihilation, (DjG)i = Gji => (DjDkG)i = Gkji etc. 
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We would like to know which actions lead to S l that are linear combinations of iterated com- 
mutators of Dj 's, i.e. Lie elements. Linearity of the passage from S(A) to 5*, implies it 
is sufficient to work with actions that are homogeneous polynomials of degree n for each 
n = 2, 3, 4 . . . separately. One difficulty is that actions are usually presented, for example, 
in the form S(A) = tr C l ^ k Ai[Aj, A k ] where C^ k are not cyclic. It takes some relabeling to 
transform to S(A) = tr S l i k Aij k with cyclic coupling tensors S^ k , in terms of which the SD 
operators are expressed in (|39p . The other difficulty is to identify those actions for which (|39p 
can be rewritten as a linear combination of iterated commutators. In what follows we carry out 
this program in part and identify a class of actions that lead to S l which are Lie elements. 

Quadratic: The most general quadratic action is 82(A) = tr C lJ AiAj. If we write it as 
5*2 = \ tr (C y + C^)Aij , then the coupling tensor 5^ = + C Ji ) is cyclically symmetric. 

By flMD, S l = 2S i W j = (C** + C^)Dj . In other words S i = (C^ + cyclic)Dj for all i are Lie 
elements, i.e. right-invariant vector fields on the group spec(Sh\) . Formally, this is true for 
arbitrary tensors C l i , though we should impose the reality condition S l i = S^* and request 
that S iJ be a positive matrix to ensure that the matrix integrals converge and lead to correlators 
satisfying G*j = Gj. In the sequel, we will work formally and not mention these reality and 
positivity conditions. 

Cubic: Motivated by the Chern-Simons example, we consider the class of cubic matrix 
models whose action is S 3 (A) = tr C ljk Ai[Aj, A k ] = tr C ljk (A ijk - A ikj ) for arbitrary C ijk . 
Writing this in the form 5 = tr (C % i k — C tk ^)Aij k = tr S l ^ k Aij k ensures that the coupling 
tensor S ijk = \(C^ + cyclic) is cyclically symmetric^. We get S l = 3S ijk D k Dj = (C ijk + 
cyclic) [D k , Dj] which are Lie elements. Notice that S l is more easily expressed in terms of the 
original C l ^ k than in terms of the coupling tensors S 1 ^ , which appear as unwelcome middlemen. 

Quartic: What is the appropriate generalization to higher degree polynomial actions, such 
that S l remains a Lie element? In [9] we showed that the Yang-Mills type of quartic ac- 
tion Sym = tr ^g^^AjWA^Ai] leads to the SD operators S l YM = Ag ik g jl [Dj, [D k , Di]} 
which are Lie elements. However, explicit calculation indicates that the more general S(A) = 
tr B l i kl [Ai, Aj][A k , Ai] leads to SD operators which are not Lie elements for some . On 
the other hand, a quartic generalization of the CS action is S4 = tr C % ^ kl Ai[A^ [A k , Ai\] . 
Furthermore, the YM action is a special case of this: using cyclicity of the trace, 

tr g ik g jl [A i , A,][A k , A t ] = tr g ik g jl (A im - A i[kl]j ) = tr g lk g jl A t [A v [A k , A t ]]. (40) 

Thus S4 reduces to Sym if C % i kl = g lk g^ 1 ■ This motivates us to check whether the SD operators 
corresponding to S4 are Lie elements for arbitrary . Write S4 as 

5 4 (A) = tr C^ k \A ijkl - A ijlk - A iklj + A ilkj ) = tr (C^ - C il ^)A m (41) 

and define S ijkl = C ij ^ - C il ^ . Then 5 4 = tr S ijkl A ijkl where the coupling tensor S ijkl = 
(§ijkl _|_ C y C ii c ) i s cyclically symmetric. Using (|39p we read off the SD operators S % = (S^ kl + 
cyclic)Di k j . After a lot of relabeling and simplification they can be written as Lie elements 
S i = \pijkl + cyclic) [[^^ D k },Dj}\ Thus in a sense, 5 4 = tr C^A^Aj, [A k ,A t ]] is the proper 
generalization of Sym while preserving the property that S l be right-invariant vector fields. 

Quintic: We begin to see a pattern to a class of actions that lead to SD operators that are 
Lie elements. 5s = tr C^ klm Ai[Aj, [A k , [Ai,A m ]]] is the obvious quintic candidate. After some 

10 The notation C l[]k] = C l]k - C ikj . 
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relabeling, we find 

S 5 (A) = tr S ijklm A ijklm , where S ijklm = C ijk[lm] - C ijm[kl] - C imm + C im,b '* ] (42) 

so that the cyclic coupling tensor is S ijklm = l(S ljklm + cyclic) and S 5 = tr S ljklm A ijk i m . 
SI = (S l i klm + cyclic)D m ;fcj after some simplification become 

S i = {c vkim + cyclic) [[ [£ , m) D^D^Dj], (43) 

which are Lie elements. While the relation of C^ klm to cyclic coupling tensors m is non- 
trivial, the SD operators are simply expressed in terms of C l ^ klm . 

Sixth degree: For the sixth degree action Sq(A) = tr C^ klmn Ai\Aj, [A k , [A h [A m ,A n ]}}}, 
the cyclic coupling tensor turns out to be S^ klmn = l(S l i klmn + cyclic) where 

gijklmn £tijkl[mn\ £jijkn[lm\ £~tijnk[lm] ^jinjk[lm] 

I £tijnm\kl] _|_ ^jinjm[kl] _|_ ^jinmj[kl] ^jinml[jk] (44) 

The corresponding S { = (C^ Mmn + cyclic) [[[[D n , D m ] , D t ], D k ] , Dj] are again Lie elements. 
Conjecture: Based on these examples, we conjecture that the n th degree polynomial action 

S n (A) = tr C^A h [A i2 , [A i3 , [• • • [A in ^,A ln ] • • •]]] (45) 

has SD operators that are the Lie elements 

5 n = {C h-in + cyclic)[[ . . . [[[A„, A„-J, A„_ 2 ], A„_ 3 ], • • •], A 2 ]- (46) 

We have exhibited a large class of matrix model actions (generalizing the gaussian, Chern-Simons 
and Yang-Mills ones) whose SD operators are right invariant vector fields on G. But it may 
not be an exhaustive list. We have not found the most general Lie elements 5* that arise as 
SD operators of some action, nor the class of all such actions. Finally, this process must be 
generalized to other groups; we hope to return to these questions later. 



6 Schwinger-Dyson operators of Yang-Mills theory 

Now we illustrate the above framework with the example of Yang-Mills theory. Instead of A 
matrices At, we now have the gluon field, one matrix A^(x) for each space-time point x and 
fj, = 1, ■ •■ ,d where d is the space-time dimension. The sources are replaced by ^(x). We 
obtain the factorized Schwinger-Dyson equations in terms of gluon correlations and write the 
SD operators S^(x) (j53[) as linear combinations of iterated commutators of left annihilation 
Du(x) with constant coefficient^^]. It follows (from sec. I4.2[) that the SD operators of Yang- 
Mills theory are right-invariant derivations of the shuffle-deconcatenation Hopf algebra generated 
by the sources where x and \x run over all space-time points and polarization indices. 

However, our formulation is far from complete. The fSDE obtained here must be supplemented 
by gauge-fixing and ghost contributions for a proper treatment of gauge invariance, before we can 
look for physical solutions. In [7] we have indicated how to incorporate gauge fixing and ghost 

11 Constant coefficients must be independent of the sources £ M (a;), but could depend on x and fj, which now 
play the role of the indices i, j, k of matrix models. Coefficients will be differential operators. 
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terms in the context of matrix models. Mathematically, this means the shume-deconcatenation 
Hopf algebra generated by needs to be modified. This could be done either by adding 

generators corresponding to ghosts or by passing to the quotient by an ideal as in Chen's work 
[5] (see also Tavares [6]). Only then can we recover the group of (generalized) loops via the 
spectrum of the Hopf algebra. We hope to return to these issues in future work, but restrict 
ourselves here to the pure Yang-Mills action 

S= tr J d*x^d ll A u (&>A v - d v A») - igd ll A v [A' t ,A u ] - ^- [A^, A v \ L4 M , A v \ j . (47) 
To get the SD equations, make the change of integration variable 

A^x) -► A'^x) = A^x) + / v^-^{x;x l ---x n )A fll (x 1 )---A fln (x n )dxi---dx n (48) 



in the Euclidean functional integral Z = J dAe~ NS . v are infinitesimal tensors and we work to 
linear order in them. This is not a local change of variable, but it is not disallowed by any law 
of physics. If 5S is the change in action, the Schwinger-Dyson equations relate it to the change 
in the measure in the sense of expectation values 

,5S. .If, (dA'\ V 

As in sec. [H in the limit N — » oo, the factorized SD equations can be written as S^(x)G(£) = 
G(£)£^(a;)Cr(£) . The generating series of gluon correlations is an element of the shuffle algebra 
generated by £^(x) ([dx] = dx\ ■ ■ ■ dx n where dxi is the volume element of space-time) 



g(o = J[dx}G fll ...^(x 1 ,---,x n )e i (x 1 )---e n (x n ) 



tr 

where G fll ... fln (xi, ■ ■ ■ , x n ) = lira (—A^ixt) ■ ■ ■ A^ n (x n )). (50) 

N^oo iv 

For the quadratic term on the rhs of the fSDE, we need the jacobianE! J = (det 

J = 1 + N 2 J dx[dx]d yi ■ ■ ■ dy k ■ ■ ■ dy n vp->*{x; x x --- x n )ffi::%- im * v " Un 

(xi -yi)-" (xjfc-i - yk-i)(x k - x)(x k+1 - y k+ i) ■■■(x n - y n ) 
tr tr 

{jjA vl ( yi ) ■ ■■A Uk _ 1 (y k - 1 )}{— A„ k+1 (y k+1 ) • ■■A Vn (Vn))- 
J -I f n 

— r = j [dx]Y J ^ 1 k " fMn (x k ;xi ■ ■ ■ x n )G llx ... llh _ 1 (xi ■ ■ ■ x k -i)G IXk+1 ... fln (x k+ i---x n ). (51) 
J k=i 

The infinitesimal change in action (^) is also linear in the arbitrary tensors v. Equating the 
coefficients of common tensors v leads to the fSDE 5 M (x)G(£) = G(£)£^(x)G(£) where the 
product on the rhs is concatenation. Let us define left annihilation D fl (x) (distinct from the 
covariant derivative) by its action on correlations: 

{D il {x)G) jJil ...^ n {xx, ■ ■ ■ ,x n ) = Gf iftl ... lin (x,xi, ■ ■ ■ ,x n ). (52) 



2 Hatted variables (e.g. dy^ ) are not integrated. 
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Then the Schwinger-Dyson operators are (square brackets denote anti-symmetrization) 

S»(x) = d v d^D^+ig{d v [D^D v ] + [0KcH, A,]} - g 2 [D\ [D»,D U ]}. (53) 

For example, let us show how S fl (x) for the 3-gluon term S3 = —ig tr f dxd^A^A^ , A u ] is 
obtained. The change in the cubic term is ( [dx] stands for dxdxi ■ ■ ■ dx n ) 



5S = -ig ti J[dx]d^v^{x ] x 1 ---x n )A^{x l )---A^n))[A^(x),A%x^ 
-ig tr J [dx^A^x))^...^; x x --- x n )A^ (x x ) ■ ■ ■ A^ (x n ), A v {x)] 
-ig tr J [dx}{d^A v {x))[A^x),v v ^ n {x- x x ■ ■ ■ x n )A^ ( Xl ) ■ ■ ■ A^{x n )\. (54) 
Integrating by parts in the first term we isolate the same v factor in all terms 

5S = ig tr J[dx]v» 1 ...^x-x x ---x n )A^(x x )---A^(x n )d p [A»(x),A u (x)] 

-ig tr J [dxy^Jx; x x ■ ■ ■ x n ){d v A^x))[A^ (x x ) ■ ■ ■ #"(x„),#W] 

-ig tr / [dx]v» 1 ...^(x;x 1 ■ ■ ■ x n )(d fl A u (x))[A^(x),A^(x 1 ) ■ ■ -A^(x n )]. (55) 



We are interested in (jt)- Before taking expectation values, we should pull d x and other 
coupling tensors to the left, while adding additional variables to ensure that derivatives act only 
on the appropriate fields. This ensures everything can be expressed in terms of gluon correlations 



5S = igg pu tr J [dx]v^ 1 ..., ln (x; x x ■ ■ ■ x n )^A^ (x x ) ■ ■ ■ A^(x n )[A"(x), A"(x)] 

-igg pp tr J [dxdy]v^...^(x;x x ■ ■ ■ x n ){x - y)d£ [a? (x)[A^ (x x ) ■ ■ ■ A^ (x n ), A"(y)] 
-igg pv tr / [dxdy]v^ n (x; x x --- x n )(x - y)d*J A» (x)[A» (y) , A* (x x ) ■ ■ ■ A^{ Xn % 



We take expectation values and use cyclicity to move fj, x ■ ■ ■ fi n to the right. This facilitates 
reexpression in terms of left annihilation. Moreover, we omit the common factor v and the 
integration over x,x x ■ ■ ■ x n , since this facilitates identifying the SD operators 

<^> « igg P „d*G^-^([x,x],x 1 ---x n )+igg pp J dy{x - y)d^G^-^([x,y],x x ■ ■ ■ x n ) 



+igg pv J dy(x - y)^G^-^([y, x],x x --- x n ) 
= igg pv dl {[D* 3 , D fi ]G)^ 1 '" tMn (x x ■■■x n ) 

+igg P »j dy{x-y)dl ([D»(y), DP{x)]Gf^{x x ■ ■ ■ x n ) 

+igg pu J dy[x-y)d? ([D"(x), D^y)]G)^" ^{x x ■ ■ ■ x n ). (57) 
We read off the SD operators of the 3-gluon terms in the Yang-Mills action 

S${x) = ig\d v [D^ D v ] + [d v D» - &>D", A,]} = igU,[D», D v ] + [d^D»\ A]}- (58) 
Similarly, the SD operators of the quadratic and 4 gluon terms in the Yang-Mills action (|47|) . 
S 2 = \ tr J d*xd tl A v (d tl A u - 8 U A»), 5 4 = -^tr J d*x[A p , A„][A^ A u ] 
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are S%(x) = d v {d»D v - d v D») = d y d^D^ 

and S%(x) = -g 2 g^g^[D u , [D p , D a )) = -g 2 [D\ [D»,D V ]\. 



(59) 



S4 contains no derivatives that need to be treated with care when identifying coupling tensors. 
So S%(x) can be read off from the commutator-squared Yang-Mills matrix model of eqn. (|37p . 
This completes the proof that the SD operators of YM theory are right invariant derivations of 
the sh-deconc Hopf algebra generated by £ M (x). This result is true independent of space-time 
dimension. 
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A Lyndon words and right standard factorization 

The SD operators of many interesting matrix models are elements of the free Lie algebra. We 
identified the FLA with the Lie algebra of right invariant vector fields on the group G = spec^ • 
To work with these vector fields it is useful to have a basis. Lyndon words introduced by Lyndon 
|26[ \27\ [28| [29] are interesting since they label a basis for the FLA. There are other bases for 
the FLA, such as the Hall basis [30]. To be self-contained, we summarize some facts (without 
proofs) about the Lyndon word basis for the FLA. 

Suppose we are given the alphabet • • • ,£ A with the order < £ 2 < • • • < £ A . There 
does not seem to be any physically preferred choice for an ordering of the letters. We extend 
the order on letters to the alphabetical or lexicographic order on all words in the alphabet. For 
example £ x £ 2 < t, 2 ^ 1 and ^S, 1 ^ 1 < C 2 *; 1 ^ 3 - If £ 7 < £ J we sa>@ £ 7 precedes £ J . A Lyndon 
word is one which is strictly minimal among its conjugates. Conjugates are words related by 
cyclic permutations; the conjugates of the Lyndon word t; 1 ^ 2 ^ 3 are ^ 2 ^ 3 ^ and S, 3 ^, 1 ^, 2 . It follows 
that Lyndon words £ L must be primitive, i.e. cannot be written as (£ M ) n for some word £ M 
and n > 2. In particular, a Lyndon word must be aperiodic, since otherwise it would equal 
one of its non-trivial conjugates. Equivalently, a word is Lyndon iff it precedes every non-empty 
proper right factor. That is, £ L is Lyndon iff for any factorization £ L = £ M £ with £ M and 
£ N non-empty, we have £ < £ . Letters are automatically Lyndon words. Lyndon words of 
length two are with < . There is also a recursive characterization of Lyndon words: 
£ is Lyndon iff there exist Lyndon words £ M and £ N such that £ < £ N and ^ L = £ M £ N ■ 
Of course, there may be more than one choice of M, N that do the job. The number of Lyndon 
words of length n over an alphabet of cardinality A is given by Witt's formula 



13 If I = 1112 ■ ■ ■ i n we abbreviate C X C 2 ' ' ' C" as £ J • F° r brevity we will sometimes talk of the word L when 
we mean the word £ L . For instance I < J really means £ 7 < £ J . 



l(n,A) 



1 



V]/x(d)A n / d , where the Mobius function 

d\n 

0, if d has a repeated prime factor; 

1, if d = 1; 

(— l) fc , if d is a product of k distinct primes. 



n 




(60) 
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For example, the numbers of Lyndon words of lengths 1,2 and 3 are given by 1(1, A) = 
A, 1(2, A) = i(A 2 - A), 1(3, A) = ±(A 3 - A). 

The right standard factorization of a Lyndon word L, \L\ > 1 is the unique factorization 
L = MN where M and N are Lyndon words such that N is of maximal length. In particular, 
M must be of length at least one. It follows that M < L = MN < N . For clarity we will denote 
the right standard factorization by L = M-N . For example, £ 1 £ 2 = £ 1 -£ 2 and C 1 ^ 3 ? 3 = C 1 ? 3- ? 3 - 

A.l Lyndon word basis for free Lie algebra 

We will specify a basis labeled by Lyndon words, for the free Lie algebra (FLA) generated by left 
annihilation D^, 1 < i < A. The FLA consists of linear combinations of iterated commutators 
of Di . Here we are implicitly thinking of the FLA as embedded in the free associative algebra 
generated by Di . Elements of the FLA are called Lie elements or Lie polynomials. A Lie element 
has a definite degree d if it is a homogeneous polynomial of degree d in the free associative 
algebra. Finding a basis for the FLA is complicated because anti-symmetry and the Jacobi 
identity relate many different Lie elements. This problem was solved |27[ [20] using the right 
standard factorization of Lyndon words. The basis elements will be called Di L \ where ^ L runs 
over all Lyndon words. Dnj) is a certain iterated commutator of Di 's where are the letters 
of the Lyndon word £ L . We use the notation -D(l) to distinguish this iterated commutator from 
the word Dl = D\ x ...i n = D^ ■ ■ ■ D[ n which is not a Lie element for n > 1 . 

Letters are Lyndon words, and the corresponding degree-one basis elements are Di , 
which are independent by definition. Given a Lyndon word L with right standard factorization 
L = M ■ N , we associate the recursively defined Lie element, Dr L \ = [-D(m)> -O(tv)] to the Lyndon 
word L. To express -D(l) as an iterated commutator of left annihilation operators, we need to 
apply this rule recursively till M and N are both single letters. The degree of the basis element 
-D(l) is equal to the length \L\. -D(l) is not defined when L is not a Lyndon word, though it is 
sometimes convenient to define -D(l) = if L is not Lyndon. 

The Lyndon words of length two are < j- Using the right standard factorization 

CC J = C " £ 3 f° r i < i, we get the Lyndon basis -D(y) = [Di,Dj] for i < j. There are clearly 
tj(A 2 — A) of these basis elements. The restriction i < j is explained by the antisymmetry of 
the commutator. 

For an alphabet of two letters (A = 2), there are two Lyndon words of length three. Their 
right standard factorizations are 

ew=^-^ 2 with e < eee < ee, 
and eee=ee-e with ee < eee < e- m 

The corresponding Lyndon basis elements are D(n2) = [D\, [Di,^]] and -D(i22) = [[Di, D2], D2] ■ 
We see that after accounting for anti-symmetry and the Jacobi identity, there are only two inde- 
pendent Lie elements of degree three for a two letter alphabet. For an alphabet of length A = 3, 
there are 8 Lyndon words of length three. Six of them involve only two of the letters each and 
can be obtained from the previous example. We list the 8 Lyndon basis elements 

[D 1 ,[D 1 ,D 2 ]}; [[£>i,D 2 ],D 2 ]; [D 1} [D 1} D 3 ]]; [[D u D 3 ], D 3 ]; 
[D 2 ,[D 2 ,D 3 ]}; [[D 2 ,D 3 ],D 3 ]; [D X ,{D 2 ,D 3 }}; [[D 1} D 3 ], D 2 }. (62) 
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The right standard factorization of the corresponding Lyndon words may be read off, for example 

£132 = £13 . £2 and £123 = £l . £23 

The structure constants of the FLA in the Lyndon basis are defined as [Drr,\, Dr^\] = 
c lmD(n)i where L,M,N are Lyndon words and \L\ + \M\ = \N\. We know of no simple 
general formula for M , but give some examples. 

D(ij), if * < j; 
[Di,Dj] = { D (ji) , if j < i; 

0, if % = j. 

1, if K = ij with i < j; 
Thus for \K\ = 2, eg = <j -1, if K = ji with j < i; (63) 

otherwise. 

Another example is (for j < k so that jk is Lyndon) 

{D(ijk), if i<jk; 

D (jki) ~ D (kji), if jk < i and k < i and k < i; 

D(j ki ) + D(j ik ), if jk < i and k < i and ji < k; ^ ' 

—D(jf.{\, if jk < i and i < k. 

Another simple example: [-D(m)>-0(jv)] = Dr MN \ if MN = M ■ N is the right standard factor- 
ization of the Lyndon word MN . For this to be the case, it is necessary that M < MN < N . 
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